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We employed in-situpu\sed laser deposition (PLD) and angle-resolved photoemission spectroscopy 
(ARPES) to investigate the mechanism of the metal-insulator transition (MIT) in NdNiOs (NNO) 
thin films, grown on NdGaO3(110) and LaA 103 ( 100 ) substrates. In the metallic phase, we observe 
three dimensional hole and electron Fermi surface (FS) pockets formed from strongly renormalized 
bands with well-defined quasiparticles. Upon cooling across the MIT in NNO/NGO sample, the 
quasiparticles lose coherence via a spectral weight transfer from near the Fermi level to localized 
states forming at higher binding energies. In the case of NNO/LAO, the bands are apparently 
shifted upward with an additional holelike pocket forming at the corner of the Brillouin zone. We 
find that the renormalization effects are strongly anisotropic and are stronger in NNO/NGO than 
NNO/LAO. Our study reveals that substrate-induced strain tunes the crystal field splitting, which 
changes the FS properties, nesting conditions, and spin-fluctuation strength, and thereby controls 
the MIT via the formation of an electronic order parameter with Qaf ~(1/4, 1/4, l/4±(5). 


INTRODUCTION 

The metal-insulator transitions (MITs) in complex ox¬ 
ides are often precursors to exotic ground states such as 
high temperature superconductivity, colossal magnetore¬ 
sistance, and different types of charge-, spin- and orbital- 
ordered states m- In this context, rare earth nickelates 
(RNiOs, where R = rare earth elements) can be viewed as 
model systems, which exhibit temperature-driven MITs 
nig. Novel routes to engineering and controlling the 
MIT in thin films and heterostructures of RNiOs further 
offer the possibility to control the MIT with thickness via 
strain using different substrates EHISI- Interestingly, the 
MIT in RNiOs, which can be tuned by changing R mg, 
by applying hydrostatic pressure [HI [ig , or by epitaxial 
strain engineering in thin films ISHlolIle], revealing that 
there are diverse possibilities for controlling correlation 
strength in these systems. 

The origin of the MIT in nickelate thin films is still 
unknown, and various explanations have been proposed, 
including magnetic order, charge order, orbital order, lat¬ 
tice distortions, or their combinations [iMi]. Experi¬ 
mentally, it is reported that the MIT temperature coin¬ 
cides with the Neel temperature T/v depending on the 
ionic radius of R HE]. More recently, various groups 
have claimed that the MIT originates from the forma¬ 
tion of charge-order due to NiOe octahedral distortions 
[2g or valence fluctuations between Ni^^'^ and Ni^”'^ ions 


[26f[^ . An observed isotope effect supports this scenario 
[30] . Eurthermore, more recent transport measurements 
on NdNiOs thin films show non-Eermi-liquid behavior 
dominating in the samples which lie at the border of the 
onset of the magnetic ordering and MIT [g. The result 
is suggestive of the presence of a ‘hidden’ quantum criti¬ 
cal point as a function of strain [g . Interestingly, optical 
spectroscopies have identified that the MIT is associated 
with a spectral weight transfer from the mid-infrared re¬ 
gion to a higher energy scale HUEll. 

Recent ARPES studies of LaNiOs thin films show 
that the ES topology changes from three- to quasi two- 
dimensional with decreasing film thickness isaisg. In 
addition, an observed ES superstructure in LNO films 
indicates some instabilities in the metallic phase [32] . 
Based on these results, it is proposed that tunable charge 
order, spin-density-wave fluctuations, and dimensional¬ 
ity might be responsible for thickness dependent MIT 
in LaNiOs and PrNiOs films isaisiES]. However, the 
origin of the possible ES instability and the mass renor¬ 
malizations responsible for the insulating state in RNiOs 
are still highly debated [22H2^ [311 [36ll4Q] . In order to 
identify the microscopic origin of the MITs and the role 
of different orders, it is crucial to understand how spec¬ 
tral functions and electronic interactions evolve across 
the transition and for different substrates. Therefore, 
in this paper, we use high-resolution ARPES to directly 
map the evolution of the electronic structure across the 
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MIT in NdNiOs (NNO) thin films grown on NdGaOs 
(NGO) and LaAlOs (LAO) substrates. 


SAMPLE PREPARATION, 
CHARACTERIZATION, AND METHODS 

Sample preparation. The 20 unit cells (uc) of NNO 
thin films were grown on NGO and LAO single-crystal 
substrates by pulsed laser deposition (PLD). A sintered 
stoichiometric NdNiOs pellet was used as an ablation tar¬ 
get. A Nb-doped yttrium-aluminum-garnet (NdiYAG) 
laser was used in its frequency-quadruple mode (266 nm) 
at a repetition rate of 2 Hz and an energy of about 
90 mJ/pulse. The substrates were heated by passing 
current through silicon wafers mounted underneath. The 
temperature of each substrate (monitored by an external 
pyrometer) was set to about 730^G and the pressure of 
the oxygen was about 0.1 mbar during the deposition. 
After deposition, each sample was cooled down to room 
temperature while maintaining the same oxygen pressure 
for about one hour. 

Characterization. In Fig. 1 we summarize the struc¬ 
tural and transport properties of the samples used in 
the ARPES experiments. Fig. 1(a) shows typical refiec- 
tion high-energy electron diffraction (RHEED) oscilla¬ 
tions during the growth of NNO thin films on NGO (for 
NNO/LAO, see Fig. S4 of supplementary information), 
indicating epitaxial growth with the potential for unit¬ 
cell thickness control. The corresponding RHEED pat¬ 
tern is shown in Fig. 1(b). Based on the growth rate 
extracted from the RHEED oscillations, the films in our 
study have total thicknesses of about 20 uc (total deposi¬ 
tion time of 240 s). In Fig. 1(c), a sharp diffraction pat¬ 
tern and some surface reconstruction derived spots are 
observed in a low-energy electron diffraction (LEED) im¬ 
age taken at 211 eV electron energy, confirming the well- 
ordered, high-quality single crystalline surfaces of the de¬ 
posited thin films. The electrical resistance of the films 
was measured by the four-point van der Pauw method 
using a physical property measurement system. The re¬ 
sistance of 20 uc NNO films on NGO(110) and LAO(100) 
substrates (the actual samples on which ARPES mea¬ 
surements were performed) is shown in Fig. 1(d) as a 
function of temperature. The MIT of nominally ^ 140 K 
is about 60 K lower than in the bulk NNO samples, con¬ 
sistent with previous reports [a E]. As can be seen in 
Fig. 1(d), despite the large variation in the strain, a sim¬ 
ilar metallic behavior from 300 K down to 170 K is ob¬ 
served in samples grown on both the substrates. How¬ 
ever, for NNO/NGO the resistance begins a step-like in¬ 
crease below about 155 K that extends until about 120 K. 
The resistance continues to increase gradually even at a 
lower temperature until a sharp increase in resistance is 
observed at around 20 K. No measurable hysteresis is ob¬ 
served down to 2 K, indicating a very sharp transition to 
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FIG. 1. Characterization of epitaxially grown NNO thin films 
on NGO(llO) and LAO(IOO) substrates, (a) Typical RHEED 
intensity oscillations and (b) corresponding RHEED pattern 
of NNO/NGO. (c) LEED pattern of NNO/NGO taken at 
211 eV electron energy, (d) Electrical resistance versus sam¬ 
ple temperature during heating and cooling cycles, as well as 
the photoemission intensity of NNO/NGO (hu = 90 eV) in¬ 
tegrated near Ef (from -0.2 to 0.1 eV) plotted along cut^l 
(as the location is marked in the upper panel of Fig. 2a) on 
the right scale, (e) X-ray diffraction pattern {0 — 20 scan) of 
NNO/NGO and NNO/LAO samples taken with Cu Ka radi¬ 
ation showing the substrate peaks. The peaks from the NNO 
films are broad due to the very low thickness (only 20 uc). 


a completely insulating state in this very low tempera¬ 
ture range. On the other hand, resistance measurements 
of NNO/LAO films show a transition at around 60 K 
with broad hysteresis, which does not appear to saturate 
at the lowest temperature of 2 K. The temperature de¬ 
pendence of the resistance across the MIT, including its 
hysteretic behavior, correlates with the changes in quasi¬ 
particle spectral weight near the Fermi level seen in the 
ARPES data [Fig. 1(d)], which will be explained in more 
detail in the Results section. 

We further characterized the films using x-ray diffrac- 
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tion at room temperature with Cu Ka radiation. Repre¬ 
sentative 0 — 20 scans of both the samples are shown in 
Fig. 1(e). No impurity phases have been detected and the 
position of the Bragg peaks are consistent with previous 
reports mm- The diffraction peaks observed at about 
47^ and 48.2^ are from (002) planes of the NGO and LAO 
substrates, respectively. A broad peak centered at about 
47.5*^ originates from the NNO thin films [as marked in 
Fig. 1(e)]. This broadening of the film peaks is expected 
because our films are only 20 uc thick. The full width 
at half maximum of the film is found to be about 0.15° 
in the rocking curve measurements (not shown), confirm¬ 
ing the good quality of the NNO films (slight broadening 
may be due to the presence of dislocations [9]). By com¬ 
paring the in-plane pseudocubic lattice constants of the 
NGO (a = 3.858 A) and LAO (a = 3.794 A) substrates 
with NNO (a = 3.803 A), the NNO films are nominally 
under tensile strain of +1.4% and compressive strain of 
-0.3% when deposited on NGO and LAO substrates, re¬ 
spectively. 

ARPES measurements. The freshly grown NdNiOs 
thin film samples were studied in-situ using ARPES. 
The ARPES measurements were performed at the Sur¬ 
face/Interface Spectroscopy (SIS) X09LA beamline of the 
Swiss Light Source located at the Paul Scherrer Institute 
in Villigen, Switzerland [41] . The end station is equipped 
with a VG-Scienta R4000 electron energy analyzer and a 
six axis liquid He-cooled manipulator (GARVING). We 
used photon energies from 60 eV to 115 eV with both 
circular and linear photon polarizations. Spectra were 
acquired over a temperature range from 20 to 200 K. 
The energy and momentum resolutions were set to about 
30 meV and « 0.009/0.019 A“^ (parallel/perpendicular 
to the analyzer slit), respectively. The binding energy 
scale was calibrated with a polycrystalline copper refer¬ 
ence sample in direct electrical and thermal contact with 
the film. The base pressure of the UHV chamber during 
the measurements was below 5x10“^^ mbar. Our results 
were reproduced on several thin film samples. 

Theoretical calculations. The theoretical calcula¬ 
tion of the complex self-energy effects is performed within 
the momentum dependent density fluctuations (MRDE) 
model [42l [43]. We first derive a three-orbital tight- 
binding model to capture the essential non-interacting 
band structure. We include nickel da, 2 _^ 2 , dz ‘2 orbitals 
whose related bands are observed to cross Ep in all our 
samples, and also the dxy orbital, which crosses Ep for 
the NNO/LAO case. The substrate dependent crystal 
field effect is parametrized by tunning the onsite poten¬ 
tial difference of the dxy orbital with respect to its eg 
counterparts, as well as by changing the inter-orbital hop¬ 
ping between them. This causes not only that the dxy 
orbital shift toward to Ep, but also it helps to change 
the ES topology of the other two Cg orbitals, in accor¬ 
dance with experiment. We note that the Ni dxy orbital 
is strongly hybridized with O p-orbitals, and thus we fit 


to the overall bandwidth. The tight-binding fitting to 
the density-functional theory band structure is given in 
the Supplementary Information. 

The electronic interactions in transition metal oxides 
are dominated by density-density correlations which in¬ 
clude spin, charge and orbital fluctuations. The den¬ 
sity fluctuation susceptibilities in the particle-hole chan¬ 
nel represent the joint density of states (JDOS), which 
is calculated by convoluting the corresponding Green’s 
function over the entire Brillouin zone (BZ). The many- 
body corrections are captured within the random-phase 
approximation (RPA) by including the onsite multiband 
components such as intra-, inter-orbital Hubbard inter¬ 
actions (t/, V), Hund’s coupling Jh^ and the pair-change 
term J' . Their values are found to be same for both sam¬ 
ples with U=2 eV for (ia,2_^2,1.5 eV for the d^'^ orbital 
and 0.8 eV for the dxy orbital. The inter-orbital inter¬ 
action is V =1.5 eV, Hund’s coupling J^=0.8 eV, and a 
pair-hoping term J' =0.33 eV for all orbitals. 

We calculate the electronic self-energy due to the cou¬ 
pling of the density-fluctuations to the electronic states 
within the MRDE method [421444] . 
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Here k and uj are the quasiparticle momentum and fre¬ 
quency, and q and uOp are the bosonic excitation mo¬ 
mentum and frequency, respectively, m, n, s, and t are 
the orbital indices, and v is the quasiparticle band index. 
The index i stands for different types of fluctuatuion such 
as spin, charge, etc. (Ibz is the electronic phase space 
volume. and Up are the fermion and boson occupa¬ 
tion numbers, and is the corresponding tight-binding 
bands. The vertex correction encodes both 

the angular and dynamical parts of the vertex, which 
are combined to obtain r^„ „(k,q) = </>k-q,sK-q,td - 
9IlAni(k — q,cd)/du;)o, where 0 is the eigenvector. 
Knh i^ fluctuation-exchange potential gen¬ 

erated by the exchange of spin, valence (charge), and 
orbital between different electrons: V^l^^(q,cjp) = 

{irii/2)\UiXi{q^^p)Uj\^^^^ where i stands for spin and 
charge (and orbital) components with r]i = 3, 1, respec¬ 
tively. Ui are the Goulomb interaction matrices, defined 
in the orbital basis, whose components are 7/, V, Jh, and 
J' as defined in the supplementary material. Xi are the 
RPA spin and charge susceptibilities. Finally, the self¬ 
energy dressed Green’s function becomes G~^ = — 

where E is the total self-energy matrix. We compute all 
self-energies self-consistently, until the total self-energy 
converges, as described in Refs. [421444] . 
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FIG. 2. Fermi surfaces (FSs) measured at different photon 
energies. (Upper panels) The FS maps of 20 uc NdNiOs thin 
films grown on NdGaOa, measured at T = 200 K with (a) 
hiy = 90 eV (k^ 0), (b) 80 eV, (c) 75 eV, and (d) 65 eV 

(k^ tt). The photoemission intensity has been integrated 
over ±5 meV. The lower panels show the calculated FSs at 
the roughly corresponding kz values. 


RESULTS 
I - NNO/NGO 


Figures 2(a)-(d) show ARPES data of NdNiOs (20 uc) 
grown on NdGaO3(110). The maps are 2D cuts through 
the Fermi surface in the k^^-k^ plane at different val¬ 
ues accessed using 90, 80, 75 and 65 eV photon energies 
(a more complete set of maps is presented in Fig. S4 of 
the Supplementary Material). The data were collected 
in the metallic phase at 200 K. Note that all the experi¬ 
mental ARPES data presented here are measured in the 
second Brillouin zone (BZ). At ku = 90 eV [Eig. 2(a)], we 
observe a nearly circular ES pocket at the center of the 
BZ (F-point), which becomes smaller when viewed with 
photon energies of 80 eV [Fig. 2(b)] and 75 eV [Fig. 2(c)]. 
At the same time, clear features appear centered around 
the zone corners. For ku = 65 eV, the photoemission 
intensity of the circular FS pocket at the zone center dis¬ 
appears completely, and a sharp band around the zone 
corner [A-point in Fig. 2(d)] forms a larger squarelike FS 
pocket. Therefore, our photon energy dependent mea¬ 
surements demonstrate the strong three-dimensional na¬ 
ture of these FS pockets. The measured FS topology is in 
good agreement with the theoretical model calculations 
(lower panel of Fig. 2) at the kz points corresponding 
with the photon energies. 

Next we show the nature of the FS pockets by studying 
the low energy band dispersions. From our detailed ku 
dependent FS maps, we focus here on 90 eV [cut 7^1 in 
3(a) and cut#2 in 3(b)] and 65 eV [cut#3 in 3(c), cut#4 


in 3(d)] photon energies, which are close to the F and Z- 
points, respectively. Figs. 3(a)-(d) show cuts along four 
different high-symmetry lines M — F — M, X — F — X, 
A — F — A, and A — R — A, respectively, which we refer to 
throughout the text as cuts7^1-4. In the metallic phase 
(175 K), the data reveal that the small pocket at the F- 
point is electron-like, while the larger pocket centred at 
the BZ diagonal (A—point) is hole-like [45f[47] . Theoret¬ 
ical calculations find that the two bands have predom¬ 
inantly dx 2 _y 2 and d^z‘^-r‘^ characters, respectively (see 
Fig. SI). 

As we cool the sample, no significant change in the 
band dispersion and intensity is visible for temperatures 
down to 125 K. However, an abrupt decrease in quasipar¬ 
ticle spectral weight is observed when the temperature 
crosses below 100 K. Warming the sample reverses the 
behavior and restores spectral weight close to the Fermi 
level [Eigs. 3(a,b)] following a hysteretic behavior that is 
also reflected in transport data [see Eig. 1(d)]. The loss of 
spectral weight near the Eermi level occurs for both the 
electron [Eigs 3(a,b)] and the hole [Eigs 3(c,d)] pockets. 
We visualize the redistribution of spectral intensity be¬ 
tween the metallic and the insulating state by subtracting 
the low temperature data from the corresponding high 
temperature ones, as plotted in Eigs. 3(e-h). The spec¬ 
tral weight transfer from the Ep region to higher BE 
is clearly visible in the entire BZ, which is found to be 
strongly momentum-dependent, unlike in a strongly cor¬ 
related Mott scenario. 

In Eig. 4 we provide more comprehensive analysis of 
the spectral weight transfer in the momentum distri¬ 
bution curves (MDGs) and energy distribution curves 
(EDCs) at different temperatures. In Eigs. 4(a,b) and 
4(e,f), the MDGs of electron and hole bands extracted 
from the panels (a,b) and (c,d) of Eig. 3 are presented. 
The corresponding EDGs at the Eermi momentum (kp) 
are shown in Eigs. 4(c,d) and 4(g,h). Note that the data 
points shown earlier in Eig. 1(d) were obtained by in¬ 
tegrating the EDGs in Eig. 4(c) from -0.2 to 0.1 eV. 
These show a step-like, hysteretic change in spectral in¬ 
tensity centred near the transport-measured transition 
temperature of about 140 K. All MDGs at 175 K exhibit 
two well-defined peaks with substantial spectral weight 
at the Eermi level. There is virtually no change in the 
MDG peak positions as the temperature is reduced; how¬ 
ever, the peak intensity is significantly suppressed below 
125 K relative to higher temperatures. In all the cuts, 
the EDGs at 175 K show a high intensity quasiparticle 
peak extending from the Eermi level up to Ep — 0.2 eV. 
There is no drastic change upon lowering the tempera¬ 
ture up to 125 K; however, a significant decrease in in¬ 
tensity is observed at around 100 K. The quasiparticle 
spectral weight continues to decrease as the temperature 
is lowered to about 75 K after which point the inten¬ 
sity loss appears to saturate (i. e., the EDGs for 75 K 
and 20 K are almost overlapping). Most intriguingly. 
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FIG. 3. ARPES spectra from 
NNO/NGO as a function of temper¬ 
ature. Spectra measured with hiy = 
90 eV (F-point) at various sample 
temperatures are plotted in the up¬ 
per two panels [cut#l (M — F — M) 
in (a) and cut#2 (X-F-X) in (b)]. 
Gorresponding spectra obtained at 
the Z-point, measured with hu = 
65 eV, at various sample tempera¬ 
tures, are plotted in the lower two 
panels [cut#3 (A — Z — A) in (c) and 
cut#4 {A — R — A) in (d)]. In (e-h), 
the difference band dispersion plots 
are obtained by subtracting 175 K 
from 75 K for all the cuts. The lo¬ 
cations of the cuts are marked in 
Figs. 1(a) and 1(d). 
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FIG. 4. MDGs and EDGs compared 
as a function of sample temperature 
during cooling and heating cycles be¬ 
tween 175 K and 20 K. (a-b) MDGs 
at Ey and (c-d) EDGs at /cp from the 
spectra in Eig. 2(a,b) {hv — 90 eV, 
F-point). (e-f) Analogous MDGs 
and (g-h) EDGs from the spectra in 
Eig. 2(c,d) (hz/ = 65 eV, Z-point). 
The insets of (c) and (d) compare the 
EDGs over a wider energy range, re¬ 
vealing the spectral weight transfer 
from near the Eermi level to higher 
binding energy. 


the intensity below 0.2 eV increases with decreasing tem¬ 
perature. Earlier angle-integrated photoemission studies 
of Ndi-a^Sm^^NiOs show similar behavior of the spectral 
weight transfer from near Ey region to higher binding en¬ 
ergy region [36]. Our observation is also consistent with 
an optical study, where the MIT is argued to be linked 
to a spectral weight transfer from the Drude region (co¬ 
herent state) to the deeper energv fTTl [31] . 

II - NNO/LAO 

In order to gain insight into the origins of the MIT in 
NNO and the accompanying spectroscopic changes seen 
by ARPES, we have leveraged the natural advantages 
of thin films as tunable systems, exploiting, in partic¬ 
ular, different substrates to change the epitaxial strain 
and thereby modify the electronic properties of NNO. 
Earlier studies have shown that substrate-induced strain 
alters the MIT, but the specific mechanism by which it 
influences the MIT has remained unknown. The temper¬ 
ature dependence of the resistivity of NNO/LAO(001), 
where the strain is compressive (~ —0.3%), is notably 


different from that of the tensile-strained (« +1.4%) 
NNO/NGO(110) films studied above. Namely, our 20- 
uc thick NNO/LAO films show metallic behavior down 
to 60 K with broad hysteresis in the MIT that is not 
saturated down to the lowest measurement temperature 
of 2 K, whereas the NNO/NGO samples have an MIT 
temperature of about 140 K and a narrower hysteresis of 
about 30 K [see Fig. 1(d)]. While the strain-dependent 
behavior of the MIT in NNO films varies somewhat 
among recent studies, the reduced transition tempera¬ 
ture we observe in NNO/LAO is in qualitative agreement 
with recent transport results [9]. Another study found 
that the temperature-driven MIT can even be completely 
suppressed under compressive strain, which appears to 
coincide with the disappearance of antiferromagnetic or¬ 
dering that otherwise is found in tensile-strained NNO 
|6] . These observations suggest a unique role of the com¬ 
pressive strain in suppressing the insulating phase, as well 
as magnetic ordering at low temperature. 

Reflecting their different electronic behavior, ARPES 
reveals marked differences between the electronic struc¬ 
tures of tensile-strained NNO/NGO and compressively- 
strained NNO/LAO. Figs 5(a-d) show FS maps mea- 
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(g) NNO/NGO, expt. (h) NNO/LAO, expt. (i) NNO/NGO, cal. (j) NNO/LAO, cal. 
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FIG. 5. (a-d) Fermi surface maps of 20 uc NdNiOs thin 

films grown on LaAlOa, measured at 200 K and 20 K using 
\w — 90 eV, and 65 eV. In order to make the extra hole¬ 
like pocket (around A-point) clearly visible, we have marked 
the position by dotted circle around A-point in the upper 
corner of (c,d). Theoretically calculated FSs are plotted in 
(e,f). (g,h) Band dispersions of NNO/NGO and NNO/LAO 
and (i, j) plots of the theoretically calculated ARPES spec¬ 
tral function (including self-energy) for the same cut#4. The 
corresponding non-interacting bands are overlaid in (i,j). 


sured with 90 eV (F-point) and 65 eV (Z-point) photon 
energies, which show the persistence of the metallic FS 
down to 20 K. The overall shape of the FS remains similar 
at high and low temperatures, and the quasiparticles- 
which are still evident at 20 K—sharpen as the tem¬ 
perature is decreased. Resistance measurements of our 
NNO/LAO also show the absence of a complete insulat¬ 
ing phase down to 2 K [Fig. 1(d)]. At 200 K, the FS of 
NNO/LAO contains an electron pocket centered at F and 
cuboidal hole pockets centered at the BZ corners. These 
resemble the Fermi surface components of NNO/NGO, 
except that the hole pockets in NNO/LAO are notice¬ 
ably larger and the electron pocket is slightly smaller. 
Strikingly, however, in NNO/LAO, extra pockets appear 
at the corners of the BZ [see around A-point at k^,, k^ « 
3, -1 and 3, 1 in Figs. 5(c,d)]. Such bands are not visible 
in NNO/NGO within the binding energy range extend¬ 
ing ^300 meV from Ep that we have studied. This ob¬ 
servation suggests that the strain not only causes band 
renormalization as one would anticipate in general, but it 
also changes the band ordering. With theoretical mod¬ 
eling we find that strain also changes the crystal field 
splitting, pushing the t 2 g orbitals upward, while pulling 
the Cg orbitals further down. This makes the electron- 
pocket at the F pocket shrink, while the hole-pocket at 
the Z point grows [as shown in Fig. 5(g,h) by extract¬ 
ing the peak positions of MDCs at each energy and in 
Figs. 5(i,j) by the theoretical calculations along cut#4 
{A — R — A)], and additionally a tiny hole pocket also 
forms at (tt, tt, tt) point as the dxy orbitals now cross the 
Fermi level in NNO/LAO system. 


DISCUSSION 

Our ARPES experiments reveal two important obser¬ 
vations related to the MIT. Firstly, the spectral weight 
transfer across the MIT from the quasiparticle states 
near Fermi level to the high-energy incoherent states is 
strongly anisotropic, indicating that the correlation effect 
is momentum dependent. Secondly, the substrate driven 
strain changes the FS areas, and thereby controls the FS 
nesting driven insulating gap. In order to understand 
the mechanism of MIT, we have combined our results 
with the theory. The momentum-dependent band renor¬ 
malization and spectral weight transfer is well captured 
within the MRDF theory, as demonstrated in Fig. 6 (for 
more detail, see the supplementary information). The 
model captures the intermediate coupling regime of the 
interaction, in which the effective interaction is of the 
order of the bandwidth [42]. In this limit, the electronic 
spectrum is often split into low-energy itinerant and high- 
energy localized states by coupling to spin, charge and 
other density fluctuations. The added benefit of the 
MRDF model is that it captures the full momentum and 
energy dependence of the self-energy, which smoothly in¬ 
terpolate the band and spectral weight between the itin¬ 
erant and localized states. 

The calculated self-energy-dressed spectral weight 
maps are presented in Figs. 6(a) and 6(b), and com¬ 
pared with the corresponding tight-binding band struc¬ 
ture (solid lines) for NNO/NGO and NNO/NGO, respec¬ 
tively. Considerable deviation of the calculated spec¬ 
tral function with respect to its non-interacting band is 
visible here. Both real (S') and imaginary (E") parts 
of the spectral function (which are related to the band 
renormalization and lifetime, respectively) are plotted 
in Fig. 6(c,d) at /c = 0,0,0 for both NNO/NGO (solid 
lines) and NNO/LAO (dashed lines), while its momen¬ 
tum dependence is presented in the supplementary infor¬ 
mation. For the NNO/NGO, we notice that the self¬ 
energies have higher intensity, and the corresponding 
slope in E' (determining the strength of mass enhance¬ 
ment at the Fermi level) is also larger than the case of 
NNO/LAO. E'{E < 0) shows a peak around 0.2-0.3 eV, 
while E" shows a peak where E' changes sign. These 
calculated results |42H44| are in good agreement with 
our ARPES experimental observation of the crossover 
that occurs at E^ — 0.2 eV between the strongly renor¬ 
malized band near E^' and the localized like states at 
higher energy. The computed band renormalization fac¬ 
tor Zk = (1 — 9E'(k, cc;)/9co’)“^q changes characteristi¬ 
cally throughout the BZ, and is found to be stronger in 
NNO/NGO than NNO/LAO, as discussed in the supple¬ 
mentary information. 

Secondly, our ARPES data demonstrated that the FS 
areas changes between NNO/NGO and NNO/LAO. In 
addition, for the NNO/LAO system, a band crosses Ep 
at the A-point, forming the extra hole like pockets [also 
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FIG. 6. Theoretical results of the spectral weight renormalization and self-energy of NNO/NGO and NNO/LAO. (a,b) Spectral 
weight maps along high-symmetry momentum cuts, computed with the inclusion of a many-body E(k, due to dynamical 
electron-electron interactions, and compared with the corresponding non-interacting band structure, plotted in black solid line, 
(c, d) Real and imaginary parts of self-energy at k = 0,0,0 for both NNO/NGO (solid lines) and NNO/LAO (dashed lines) 
systems, (e, f) Full spectrum of many-body spin-correlation function, (g) susceptibility compared in the low-energy region. 
(h,i) calculated Fermi surfaces for both the systems. 


can be seen in the calculated FSs as shown in Fig.6(i)]. 
In order to explain this finding, we model the crystal 
field splitting by parameterizing the relative onsite po¬ 
tential between the low-lying Cg and t 2 g orbitals, as well 
as the hybridization between them. These two param¬ 
eters lower the chemical potential for the eg orbitals, 
while lifting the t 2 g orbital closer to the Fermi level in 
NNO/LAO, as seen experimentally. Recently published 
x-ray absorption spectra also show a remarkably large 
and asymmetric energy shift in the O K-edge as a func¬ 
tion of strain [6], suggesting the same conclusion. H. K. 
Yoo, et al observed the strain dependent FS superstruc¬ 
tures in LaNiOs films [32], which may be due to the mod¬ 
ification in the low-lying orbitals and crystal field split¬ 
ting. It has been suggested that the existence of charge 
disproportionation in RNiOs is strongly related to the 
metal-insulator phase transition [26H29] . Our ARPES 
results on two differently strained films suggest the exis¬ 
tence of a tunable electronic localization. If so, then the 
number of itinerant charges, above the transition in the 
metallic state can be different for NNO films grown on 
NGO and LAO. Indeed, we find that strain-tuned crystal 


field splitting alters the FS nesting conditions, in turn in¬ 
fluencing the spin fluctuations that play a dominant role 
in the MIT. 

We further performed the calculations in order to cap¬ 
ture the full spectrum of the spin susceptibility for the 
strained NNO films, see Figs. 6(e,f). Notably, we find 
that the spin susceptibility dominates over the charge 
susceptibility in both systems. The RPA susceptibilites 
exhibit two dominant intensities at different energy and 
momentum scales. The low-energy collective modes, such 
as magnon with preferred wavevector indicate the pres¬ 
ence of electronic or magnetic order, while the high- 
energy dispersive modes, like paramagnons, gives the 
marginal Fermi-liquid-like state. The spin susceptibil¬ 
ity, in the low-energy region, is stronger in NNO/NGO 
and concentrated around the wavevector Qaf( 1/4, 1/4, 
l/4±(5)7r/a, while for NNO/LAO, it is spread out over a 
large momentum region, as shown in Fig. 6(e-g). The ori¬ 
gin of the stronger spin-fluctuation spectral weight at one 
dominant Qaf value in NNO/NGO [Fig. 6(g)] is related 
to the underlying band structure and nesting-driven FS 
instability, which is found to be weaker in NNO/LAO. 















This indicates that NNO/NGO is prone to a magnetic 
instability through dynamical fluctuations resulting in 
its stronger band renormalization found in the metallic 
state, which ultimately may account for its higher MIT. 

CONCLUSION 

In summary, we have presented the first momentum- 
resolved study of the electronic structure of NNO thin 
films across the MIT as a function of epitaxial strain (us¬ 
ing two different substrates). Both compressive and ten¬ 
sile strained films in the metallic phase showed 3D FS and 
share common elements: spheroidal electron like pockets 
centered at T and large cuboidal hole like pockets at the 
3D BZ corners. The corresponding bands are strongly 
renormalized in the metallic phase, in agreement with 
theoretical calculations. A transfer of the spectral weight 
from near Ey region to higher binding energy occurs in 
the entire BZ across the MIT. On the basis of these ob¬ 
servations, we have performed a self-energy calculation 
due to spin and charge fluctuation in the dynamical cor¬ 
relation limit, which shows that the intertwined FS nest¬ 
ing and the dynamical correlation effects are stronger in 
NNO/NGO than NNO/LAO. We find that the strain re¬ 
duces crystal field splitting in NNO/LAO, lowering the 
Cg states and lifting those of the t 2 g orbitals. This leads 
to the appearance of a new holelike FS at the A—point, 
which plays a crucial role in lowering the MIT. Our study 
provides insights into the physical mechanisms responsi¬ 
ble of driving the MIT in complex oxides, implying how 
strain, nesting, correlations, and magnetic fluctuations 
can be tuned, which may help the future development of 
rare earth nickelates in applications [9l [48H5Q] . 

ACKNOWLEDGMENT 

The authors thank technical staff at the SIS beam line 
for their cooperation and help during the experiments. 


* rsdhaka@physics.iitd.ac.in 
^ milan.radovic@psi.ch 
^ joel.mesot@psi.ch 

[1] M. Imada, A. Fujimori, and Y. Tokura, Metal-insulator 
transitions, Rev. Mod. Phys. 70 , 1039 (1998). 

[2] H. Y. Hwang, Y. Iwasa, M. Kawasaki, B. Keimer, N. 
Nagaosa, and Y. Tokura, Emergent phenomena at oxide 
interfaces. Nature Materials 11 , 103 (2012). 

[3] P. Zubko, S. Gariglio, M. Gabay, P. Ghosez, J. M. 
Triscone, Interface physics in complex oxide heterostruc¬ 
tures, Annu. Rev. Condens. Matter Phys. 2 , 141 (2011). 

[4] M. Medarde, Structural, magnetic and electronic prop¬ 
erties of RNiOa perovskites (R = rare earth), J. Phys.: 
Condens. Matter 9 , 1679 (1997). 


[5] G. Catalan, Progress in perovskite nickelate research. 
Phase Transitions 81 , 729 (2008). 

[6] J. Liu, M. Kargarian, M. Kareev, B. Gray, P. J. Ryan, A. 
Cruz, N. Tahir, Yi-De Chuang, J. Guo, J. M. Rondinelli, 
J. W. Freeland, G. A. Fiete, and J. Chakhalian, Het¬ 
erointerface engineered electronic and magnetic phases 
of NdNiOa thin films. Nature Communications 4, 2714 
(2013). 

[7] J. Liu, M. Kareev, D. Meyers, B. Gray, P. Ryan, J. W. 
Freeland, and J. Chakhalian, Metal-insulator transition 
and orbital reconstruction in Mott-type quantum wells 
made of NdNiOa, Phys. Rev. Lett. 109 , 107402 (2012). 

[8] R. Scherwitzl, S. Gariglio, M. Gabay, P. Zubko, M. Gib- 
ert, and J.-M. Triscone, Metal-insulator transition in 
ultrathin LaNiOs films, Phys. Rev. Lett. 106 , 246403 
( 2011 ). 

[9] R. Scherwitzl, P. Zubko, 1. G. Lezama, S. Ono, A. F. Mor- 
purgo, G. Catalan, and J. M. Triscone, Electric-field con¬ 
trol of the metal-insulator transition in ultrathin NdNiOa 
films. Advanced Materials 22 , 5517 (2010). 

[10] G. Catalan, R. M. Bowman, and J. M. Gregg, Metal- 
insulator transitions in NdNiOs thin films, Phys. Rev. B 
62 , 7892 (2000). 

[11] M. K. Stewart, D. Brownstead, J. Liu, M. Kareev, J. 
Chakhalian, and D. N. Basov, Heterostructuring and 
strain effects on the infrared optical properties of nicke¬ 
lates, Phys. Rev. B 86, 205102 (2012). 

[12] E. Benckiser et al.^ Orbital reflectometry of oxide het¬ 
erostructures. Nat. Materials 10 , 189 (2011). 

[13] M. Wu et aL, Strain and composition dependence of or¬ 
bital polarization in nickel oxide super lattices. Phys. Rev. 
B 88, 125124 (2013). 

[14] P. C. Canfield, J. D. Thompson, S. W. Cheong, and L. W. 
Rupp, Extraordinary pressure dependence of the met al¬ 
to-insulator transition in the charge-transfer compounds 
NdNiOa and PrNiOa, Phys. Rev. B 47 , 12357 (1993). 

[15] X. Obradors, L. M. Paulius, M. B. Maple, J. B. Tor¬ 
rance, A. 1. Nazzal, J. Eontcuberta, and X. Granados, 
Pressure-dependence of the metal-insulator-transition 
in the charge-transfer oxides RNiOs (R = Pr, Nd, 
Ndo.rLao.s), Phys. Rev. B 47, 12353 (1993). 

[16] A. M. Kaiser, A. X. Gray, G. Conti, J. Son, A. Greer, 
A. Perona, A. Rattanachata, A.Y. Saw, A. Bostwick, 
S. Yang, S.-H. Yang, E. M. Gullikson, J. B. Kor- 
tright, S. Stemmer, and C. S. Eadley, Suppression of 
near-Eermi level electronic states at the interface in 
a LaNiOs/SrTiOs superlattice, Phys. Rev. Lett. 107 , 
116402 (2011). 

[17] J. L. Garcfa-Muhoz, M. A. G. Aranda, J. A. Alonso, 
and M. J. Martmez-Lope, Structure and charge order in 
the antiferromagnetic band-insulating phase of NdNiOs, 
Phys. Rev. B 79 , 134432 (2009). 

[18] V. Scagnoli, U. Staub, Y. Bodenthin, M. Garc’ia- 
Eern’andez, A. M. Mulders, G. 1. Meijer, and G. Ham- 
merl. Induced noncollinear magnetic order of Nd^^ in 
NdNiOs observed by resonant soft x-ray diffraction, 
Phys. Rev. B 77 , 115138 (2008). 

[19] M. Medarde, A. Eontaine, J. L. Garcfa-Muhoz, J. 
Rodrguez-Carvajal, M. de Santis, M. Sacchi, G. Rossi, 
and P. Lacorre, RNiOs perovskites (R = Pr, Nd): Nickel 
valence and the metal-insulator transition investigated by 
x-ray absorption spectroscopy, Phys. Rev. B 46 , 14975 
(1992). 



9 


[20] D. Kumar, K. P. Rajeev, J. A. Alonso, and M. J. 
Martmez-Lope, Spin-canted magnetism and decoupling 
of charge and spin ordering in NdNiOs, Phys. Rev. B 
88, 014410 (2013). 

[21] V. Scagnoli, U. Staub, A. M. Mulders, M. Janousch, 

G. I. Meijer, G. Hammer 1, J. M. Tonnerre, and N. 
Stojic, Role of magnetic and orbital ordering at the 
metal-insulator transition in NdNiOs, Phys. Rev. B 73, 
100409(R) (2006). 

[22] M. J. Han, X. Wang, G. A. Marianetti, and A. J. Millis, 
Dynamical mean-field theory of nickelate super lattices, 
Phys. Rev. Lett. 107, 206804 (2011). 

[23] H. Park, A. J. Millis, and G. A. Marianetti, Site-selective 
Mott transition in rare-earth-element nickelates, Phys. 
Rev. Lett. 109, 156402 (2012). 

[24] B. Lau, and A. J. Millis, Theory of the magnetic and 
metal-insulator transitions in RNiOs bulk and layered 
structures, Phys. Rev. Lett. 110, 126404 (2013). 

[25] S. Johnston, A. Mukherjee, 1. Elfimov, M. Berciu, and G. 
A. Sawatzky, Gharge disproportionation without charge 
transfer in the rare-earth-element nickelates as a possible 
mechanism for the metal-insulator transition, Phys. Rev. 
Lett. 112 , 106404 (2014). 

[26] U. Staub, G. 1. Meijer, F. Fauth, R. Allenspach, J. G. 
Bednorz, J. Karpinski, S. M. Kazakov, L. Paolasini, and 
F. d’Acapito, Direct observation of charge order in an epi¬ 
taxial NdNiOs film, Phys. Rev. Lett. 88, 126402 (2002). 

[27] V. Scagnoli, U. Staub, M. Janousch, A. M. Mulders, M. 
Shi, G. 1. Meijer, S. Rosenkranz, S. B. Wilkins, L. Pao¬ 
lasini, J. Karpinski, S. M. Kazakov, and S. W. Lovesey, 
Gharge disproportionation and search for orbital order¬ 
ing in NdNiOs by use of resonant x-ray diffraction, Phys. 
Rev. B 72, 155111 (2005). 

[28] M. Medarde, G. Dallera, M. Grioni, B. Delley, F. Vernay, 
J. Mesot, M. Sikora, J. A. Alonso, and M. J. Martinez- 
Lope, Gharge disproportionation in RNiOs perovskites 
(R = rare earth) from high-resolution x-ray absorption 
spectroscopy, Phys. Rev. B 80, 245105 (2009). 

[29] K. Bilewska, E. Wolna, M. Edely, P. Ruello, and J. Szade, 
Evidence of charge disproportionation on the nickel sub¬ 
lattice in EuNiOs thin films: X-ray photoemission stud¬ 
ies, Phys. Rev. B 82, 165105 (2010). 

[30] M. Medarde, P. Lacorre, K. Gonder, F. Fauth, and 
A. Furrer, Giant ^®0-^®0 isotope effect on the metal- 
insulator transition of RNiOs perovskites ( R = Rare 
Earth), Phys. Rev. Lett. 80, 2397 (1998). 

[31] M. K. Stewart, J. Liu, M. Kareev, J. Ghakhalian, and D. 
N. Basov, Mott physics near the insulator-to-metal tran¬ 
sition in NdNiOs, Phys. Rev. Lett. 107, 176401 (2011). 

[32] H. K. Yoo, et al., Dimensional crossover of the elec¬ 
tronic structure in LaNiOs ultrathin films: Orbital re¬ 
construction, Fermi surface nesting, and the origin of 
the metal-insulator transition, arXiv: 1309.0710; Latent 
instabilities in metallic LaNiOs films by strain control of 
Fermi-surface topology. Scientific Reports 5, 8746 (2015). 

[33] P. D. G. King, H. 1. Wei, Y. F. Nie, M. Uchida, G. Adamo, 
S. Zhu, X. He, 1. Bozovic, D. G. Schlom and K. M. Shen, 
Atomic-scale control of competing electronic phases in 
ultrathin LaNiOs, Nature Nanotechnology 9, 443 (2014). 

[34] M. Hepting, M. Minola, A. Frano, G. Gristiani, G. 
Logvenov, E. Schierle, M. Wu, M. Bluschke, E. Weschke, 

H. -U. Habermeier, E. Benckiser, M. Le Tacon, and B. 
Keimer, Tunable charge and spin order in PrNiOs thin 
films and superlattices, Phys. Rev. Lett. 113, 227206 


(2014). 

[35] A. V. Boris et al., Dimensionality Gontrol of Electronic 
Phase Transitions in Nickel-Oxide Superlattices. Science 
332, 937 (2011). 

[36] K. Okazaki, T. Mizokawa, A. Fujimori, E. V. Sampathku- 
maran, M. J. Martinez-Lope, and J. A. Alonso, Grossover 
in the nature of the metallic phases in the perovskite-type 
RNiOs, Phys. Rev. B 67, 073101 (2003). 

[37] E. F. Schwier, R. Scherwitzl, Z. Vydrova, M. Garcia- 
Fernandez, M. Gibert, P. Zubko, M. G. Gamier, J.-M. 
Triscone, and P. Aebi, Unusual temperature dependence 
of the spectral weight near the Fermi level of NdNiOs 
thin films, Phys. Rev. B 86, 195147 (2012). 

[38] 1. Vobornik, L. Perfetti, M. Zacchigna, M. Grioni, 
G. Margaritondo, J. Mesot, M. Medarde, and P. La¬ 
corre, Electronic-structure evolution through the metal- 
insulator transition in RNiOs, Phys. Rev. B 60, 8426(R) 
(1999). 

[39] S. R. Barman, A. Ghainani, and D. D. Sarma, Govalency- 
driven unusual metal-insulator transition in nickelates, 
Phys. Rev. B 49, 8475 (1994). 

[40] M. Medarde, D. Purdie, M. Grioni, M. Hengsberger, 
Y. Baer, and P. Lacorre, A photoemission spectroscopy 
study of PrNiOs through the metal-insulator transition, 
Europhy. Lett., 37, 483 (1997). 

[41] Y. Sassa, M. Radovic, M. Mansson, E. Razzoli, X. Y. 
Gui, S. Pailhes, S. Guerrero, M. Shi, P. R. Willmott, 
F. Miletto Granozio, J. Mesot, M. R. Norman, and L. 
Patthey, Ortho-H band folding in YBa 2 Gus 07_5 films 
revealed by angle-resolved photoemission, Phys. Rev. B 
83, 140511(R) (2011). 

[42] T. Das, R. S. Markiewicz, and A. Bansil, Intermediate 
coupling model of the cuprates, Adv. Phys. 63, 151-266 
(2014). 

[43] T. Das, J. -X. Zhu, and M. J. Graf, Spin-fluctuations and 
the peak-dip-hump feature in the photoemission spec¬ 
trum of actinides, Phys. Rev. Lett. 108, 017001 (2012). 

[44] T. Das, and K. Dolui, Superconducting dome in MoS2 
and TiSe2 generated by quasiparticle-phonon coupling, 
Phys. Rev. B 91, 094510 (2015). 

[45] P. Hansmann, X. Yang, A. Toschi, G. Khaliullin, O. K. 
Andersen, and K. Held, Turning a nickelate Fermi surface 
into a cupratelike one through heterostructuring, Phys. 
Rev. Lett. 103, 016401 (2009). 

[46] R. Eguchi, A. Ghainani, M. Taguchi, M. Matsunami, Y. 
Ishida, K. Horiba, Y. Senba, H. Ohashi, and S. Shin, 
Fermi surfaces, electron-hole asymmetry, and correlation 
kink in a three-dimensional Fermi liquid LaNiOs, Phys. 
Rev. B 79, 115122 (2009). 

[47] M. Uchida, K. Ishizaka, P. Hansmann, Y. Kaneko, Y. 
Ishida, X. Yang, R. Kumai, A. Toschi, Y. Onose, R. 
Arita, K. Held, O. K. Andersen, S. Shin, and Y. Tokura, 
Pseudogap of metallic layered nickelate R 2 -ccSrajNi 04 (R 
= Nd, Eu) crystals measured using angle-resolved pho¬ 
toemission spectroscopy, Phys. Rev. Lett. 106, 027001 
( 2011 ). 

[48] J. Mannhart, and D. G. Schlom, Oxide interfaces-An 
opportunity for electronics. Science 327, 1607 (2010). 

[49] Z. Yang, G. Ko, and S. Ramanathan, Oxide electronics 
utilizing ultrafast metal-insulator transitions, Annu. Rev. 
Gondens. Matter Phys. 41, 337 (2011). 

[50] J. Shi, S. D. Ha, Y. Zhou, F. Schools, and S. Ra¬ 
manathan, A correlated nickelate synaptic transistor. 
Nature Gommunications 4, 2676 (2013). 


